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Abstract 

The observation of ultrahigh energy neutrinos (UHEi/s) has be- 
come a priority in experimental astroparticle physics. UHEz/s can be 
detected with a variety of techniques. In particular, neutrinos can in- 
teract in the atmosphere (downward-going v) or in the Earth crust 
(Earth-skimming v) , producing air showers that can be observed with 
arrays of detectors at the ground. With the Surface Detector Array of 
the Pierre Auger Observatory we can detect these types of cascades. 
The distinguishing signature for neutrino events is the presence of very 
inclined showers produced close to the ground (i.e. after having tra- 
versed a large amount of atmosphere). In this work we review the 
procedure and criteria established to search for UHE^s in the data 
collected with the ground array of the Pierre Auger Observatory. This 
includes Earth-skimming as well as downward-going neutrinos. No 
neutrino candidates have been found, which allows us to place com- 
petitive limits to the diffuse flux of UHEi/s in the EeV range and above. 



1. Introduction 

The observation of ultrahigh energy cosmic rays (UHECR) of energy 
1 — 100 EeV (10 18 — 10 20 eV) has stimulated much experimental as well 
as theoretical activity in the field of Astroparticle Physics [TJ [2] . Although 
many mysteries remain to be solved, such as the origin of the UHECRs, their 
production mechanism and composition, we know that it is very difficult 
to produce these energetic particles without associated fluxes of ultrahigh 
energy neutrinos (UHEi/s) [3]. 

In the so-called "bottom-up" models, protons and nuclei are accelerated 
in astrophysical shocks, where pions are believed to be produced by cosmic 
ray interactions with matter or radiation at the source [3]. In the so-called 
"top-down" scenarios, protons and neutrons are produced from quark and 
gluon fragmentation, a mechanism which is known to produce much more pi- 
ons than nucleons [5j. Furthermore, protons and nuclei also produce pions in 
their unavoidable interactions responsible for the Greisen-Zatsepin-Kuzmin 
(GZK) cut-off |51 E| . The flux of UHECRs above ~ 5 x 10 19 eV is known 
to be largely suppressed with respect to that at lower energies, a feature 
seen in the UHECR spectrum [8j |9] that is compatible with the interaction 
of UHECRs with the cosmic microwave background (CMB) radiation. If 
the primaries are protons, the interaction responsible for the GZK effect is 
photopion production, and the decays of the charged pions produce UHE 
neutrinos. However, their fluxes are uncertain [3], and if the primaries are 
heavier nuclei, the UHE^ yield would be strongly suppressed [10J . 

The observation of UHE neutrinos could provide important hints to the 
origin of UHECRs |111 I12j . Unlike cosmic rays, neutrinos point directly to 
the source where they were produced, without being deflected by Galactic 



and extragalactic magnetic fields. Unlike photons they travel undisturbed 
from the sources carrying a footprint of the production model. 

High energy neutrinos can be detected with a variety of techniques [131 
114] . In particular they can be observed with arrays of detectors at ground 
level that are currently being used to measure extensive showers produced 
by cosmic rays [15]. The main challenge in this technique lies in separating 
showers initiated by neutrinos from those induced by regular cosmic rays. It 
was suggested in the 1970s that this could be done at high zenith angles [16] 
because the atmosphere slant depth provides a quite large target for neutrino 
interactions. The idea is that neutrinos, having very small cross-sections, 
can interact at any point along their trajectories, while protons, nuclei or 
photons interact shortly after entering the atmosphere. The signature for 
neutrino events is thus inclined showers that interact deep in the atmosphere. 

Inclined showers were first observed in the 1960s by several groups [TTJ, 
118] . With the Surface Detector Array (SD) of the Pierre Auger Observa- 
tory [19] we can detect inclined showers and identify neutrinos with energies 
typically above 0.1 EeV. There are two ways of performing this task: 

1. Neutrinos of all flavours can collide with nuclei in the atmosphere 
and induce an extensive air shower close to the ground [20] . In this so- 
called "downward-going" neutrino channel, both charged current (CC) 
and neutral-current (NC) interactions contribute to the neutrino event 
rate. 

2. Neutrinos of tau flavour [v T ) are expected to be most sensitively ob- 
served through the detection of showers induced by the decay products 
of an emerging r lepton, after the propagation and interaction of an 
upward-going v T inside the Earth [211 122] . This "Earth-skimming" 
channel benefits from the long range of the r lepton (~ 10 km for the 
shower energies relevant in this analysis) which sets the scale of the 
effective volume. Only charged-current interactions of v T are relevant 
in this case. 

In both the Earth-skimming and downward-going channels the showers 
can be identified and separated from cosmic ray induced showers with the SD 
of the Pierre Auger Observatory if the zenith angle is large enough, typically 
larger than ~ 65°-75°. A number of properties of the shower front, mostly 
stemming from the time distribution of the shower particles, can be used to 
distinguish neutrino-induced showers. As shown in Sec. [5j even though the 
criteria to identify neutrinos in both channels being based on similar ideas 
and variables, two different analyses were designed. The main reason for 
that concerns background reduction. The Earth-skimming neutrino search 
is restricted to a very narrow angular range where the background of nucle- 
onic showers is expected to be very small. On the other hand in the broader 



angular range of the downward-going neutrino search the background con- 
tamination is expected to be larger, and the selection criteria need to be 
more restrictive. This calls for specific algorithms and methods, capable of 
optimizing the separation of neutrino-induced showers from nucleonic ones 
as will be explained later in the paper. 

In this work we review the procedure to search for UHEi^s with the SD of 
the Auger Observatory, for both the Earth-skimming and downward-going 
channels. In Sec. [2] we give a brief overview of the SD of the Pierre Auger 
Observatory. In Sec. Owe concentrate on the general strategy to search for 
UHEi/s. Sec. H] is devoted to describe the simulations of neutrino-induced 
showers crucial to establish selection criteria and to compute the exposure to 
UHEi/s which is reported in Sec. [6j In Sec. [5] we give a detailed description 
of the neutrino selection criteria. When these criteria are applied blindly to 
the data collected at the SD no candidates are found. The resulting limits 
to the diffuse flux of UHE^s are presented in Sec. [7J Finally, in Sec. [8] we 
summarize the paper and give some prospects for future observations. 

2. The Pierre Auger Observatory 

The Pierre Auger Observatory [19] is a hybrid UHECR detector com- 
bining an array of particle detectors at ground level, and 24 fluorescence 
telescopes housed in four buildings, for redundancy and calibration. It is lo- 
cated near the town of Malargiie, in the province of Mendoza in Argentina. 
In this review we focus on the surface detector array |19[ [23] which is briefly 
described in the following. 

2.1. The Surface Detector Array 

The surface detector array [23] consists of water Cherenkov detectors in 
the form of cylinders of 3.6 m diameter and 1.2 m height, each containing 
12 tonnes of purified water. Charged particles entering the station emit 
Cherenkov light which is reflected at the walls by a diffusive Tyvek liner, 
and collected by three 9-inch photomultiplier tubes (PMT) at the top sur- 
face and in optical contact with the water. The PMT signals are sampled 
by flash analog to digital converters (FADC) with a time resolution of 25 
ns. Each station is regularly monitored and calibrated in units of vertical 
equivalent muons (VEM) corresponding to the signal produced by a muon 
traversing the tank vertically through its center [23]. In Fig. [1] we show a pic- 
ture of one of the water Cherenkov stations. The stations are autonomous, 
with all their components (PMTs, local processor, GPS receiver and radio 
system) powered by batteries coupled to solar panels. Once installed, the 
local stations work continuously without external intervention. 

The SD was completed in 2008. There are ~ 1600 water stations ar- 
ranged in a triangular grid with 1.5 km spacing between them, spanning 
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Figure 1: Top panel: one of the ~ 1600 water Cherenkov stations that 
constitute the surface detector array of the Pierre Auger Observatory (fore- 
front), and one of the four fluorescence buildings housing six of the 24 flu- 
orescence telescopes (background). Bottom panel: Layout of the SD array 
with ~ 1600 water Cherenkov stations (depicted as dots), spread over a 
surface of ~ 3000 km 2 (blue area), with a distance between stations of 1.5 
km. The four fluorescence buildings at the edges of the Observatory are also 
indicated. 11 



an almost flat surface of ~ 3000 km 2 , at an approximate altitude of 1400 
m above sea level, or equivalently an atmospheric depth -X" grounc i = 880 
g cm -2 . The layout of the SD array is sketched in the bottom panel of 
Fig.EJ 

2.2. Surface Detector trigger 

The stations transmit information by conventional radio links to the 
Central Data Acquisition System (CDAS) located in Malargiie. There are 
two types of trigger conditions. A local trigger at the level of an individual 
station (second order or T2 trigger), and a global trigger (third order or 
T3 trigger). The T2 trigger condition is the logical OR of two conditions: 
either a given threshold signal (3.2 VEM) is passed in at least one time bin 
of the FADC trace -the so-called "Threshold trigger"-, or a somewhat lower 
threshold (0.2 VEM) is passed in at least 13 bins within a 3 //s time window 
(i.e. 120 bins) -the so-called "Time-over-Threshold (ToT) trigger" . The ToT 
condition was designed to trigger on signals broad in time, characteristic of 
the early stages of the development of an extensive air shower, and is crucial 
for neutrino identification as explained below. The data acquisition system 
receives the local T2 triggers and builds a global T3 trigger requiring a 
relatively compact configuration of at least three local stations compatible 
in time, each satisfying the ToT trigger, or four triggered stations with any 
type of T2 trigger [55]. With the completed array, the global T3 trigger 
rate is about two events per minute, one third being actual shower events 
at energies above 3 x 10 eV. 

3. Generalities of UHE Neutrino search 

With the SD of the Pierre Auger Observatory we can detect and identify 
UHE neutrinos in the EeV range and above |26U27t l28]. The main challenge 
from the experimental point of view is to identify neutrino-induced showers 
in the large background of showers initiated by nucleonic cosmic rays. The 
concept for identification is relatively simple. While protons, heavier nuclei 
and even photons interact shortly after entering the atmosphere, neutrinos 
can generate showers initiated deeply into the atmosphere. When consider- 
ing vertical showers, even the ones initiated by protons or heavy nuclei have a 
considerable amount of electromagnetic component at the ground ( "young" 
shower front). However, when looking at high zenith angles {9 > 75°) the 
atmosphere is thick enough (thicker than about three vertical atmospheres) 
so that the cosmic rays interacting high in the atmosphere have shower 
fronts dominated by muons at ground ("old" shower front). A neutrino 
with 9 > 75° interacting deep will present a young shower front and, conse- 
quently, can be distinguished. 
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At the SD level, young showers induce signals spread in time over hun- 
dreds of nano-seconds in a fraction of the stations triggered by the shower, 
while old showers induce narrow signals spreading over typically tens of 
nano-seconds in practically all the stations of the event. With the 25 ns 
time resolution of the FADC of the water Cherenkov stations, the distinc- 
tion between traces induced by young and old shower fronts can be easily 
accomplished. In Fig. [2] we show an example of those two types of traces. 
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Figure 2: FADC traces of stations at 1 km from the shower core for two 
real showers of 5 EeV. Left panel: shower arriving in the early stages of 
development ("young" shower). Right panel: "old" extensive air shower 
(9 ~ 80°). 

With this simple idea, we can search for two types of neutrino-induced 
showers at the surface detector array of the Pierre Auger Observatory, 
namely: 

1. Earth-skimming showers induced by tau neutrinos {y T ) that travel in 
the upward direction with respect to the vertical to ground. v T can 
skim the Earth's crust and interact relatively close to the surface in- 
ducing a tau lepton which escapes the Earth and decays in flight in 
the atmosphere, close to the SD. 

Typically, only Earth-skimming v T -induced showers with zenith angles 
90° < 9 < 95° may be identified. 

2. Showers initiated by any neutrino flavour moving down at large angles 
with respect to the vertical at ground that interact in the atmosphere 
close to the surface detector array. We include here showers induced 
by v T interacting in the mountains surrounding the Pierre Auger Ob- 
servatory. Although this latter process is exactly equivalent to the 
"Earth-skimming" mechanism, it is included in this class because such 
showers are also going downwards. In the following we will refer to all 
these types of showers as "downward-going" i/-induced showers. 

In this review we restrict ourselves to downward-going ^-induced show- 
ers with zenith angles 75° < 9 < 90°. 
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In Fig. [3] we show a pictorial representation of the different types of 
inclined showers that can be detected. 



1) Regular proton shower. 




2) Deep Down-going u shower 



4) Down-going v T 
interacting in the 
mountains 



Muonic component of the showe 



■% 



^^3) Up-going Earth-skiming v T shower 



Figure 3: Pictorial representation of the different types of inclined show- 
ers that can be detected at the surface detector array of the Pierre Auger 
Observatory. (1) An inclined shower induced by a proton interacting high 
in the atmosphere whose electromagnetic component is absorbed and only 
the muons reach the detector. Inclined showers presenting significant elec- 
tromagnetic component at the detector level: (2) a deep downward-going 
^-induced shower; (3) an Earth-skimming v T interacting in the Earth crust 
and producing an upward-going r lepton that decays in flight and induces 
a shower in the atmosphere; and (4) a v T interacting in the mountains, pro- 
ducing a downward-going r lepton that decays close to the detector and 
initiates a shower. 



4. Simulation of neutrino showers 

Monte Carlo simulations of neutrino-induced showers are crucial to es- 
tablishing identification criteria and computing the acceptance of the SD to 
UHEi/s. The whole simulation chain is divided in three stages: 

1. High energy processes: 

a) The i/-nucleon interaction in the atmosphere for downward-going 
neutrinos is simulated with HERWIG [29] . 

The output of HERWIG includes the types, energies and mo- 
menta of the secondary particles produced for both charged (CC) 
and neutral current (NC) neutrino interactions (see Fig. U]for a 
pictorial summary of all the channels considered in this work). 

b) In the case of v T CC interactions, the r lepton propagation in the 
Earth and/or in the atmosphere is simulated with a dedicated, 
fast and flexible code which allows us to easily study the influence 
on the outgoing r lepton flux of different v T interaction cross 
sections, r energy loss models, etc. The simulation of the decay of 
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the r (when necessary) is performed with the TAUOLA package 

EDI. 



2. Shower development in the atmosphere: 

The AIRES Monte Carlo code [31] is used to propagate the particles 
produced in a high energy v interaction, or in the decay of a r lepton. 
The types, energy-momenta and times of the particles reaching the SD 
level are obtained. 

3. Surface detector array simulation: 

This is performed with the Off line software [32j . Firstly particles 
reaching a surface detector station are injected into the station, and 
with the aid of GEANT4 [33] the amount of Cherenkov light produced 
in water is calculated. Then the FADC traces of the PMT signals are 
obtained, and the total signal due to the particles entering the station, 
as well as several quantities characterizing the FADC trace which will 
be relevant for neutrino identification are computed (see below). Also 
both the local trigger condition (T2 - either Threshold or ToT), and 
the global trigger condition (T3) are applied to the simulated events 
in the same way as for collected data. 

The phase space of the simulations - namely, neutrino energy, zenith 
angle of incidence, interaction depth in the atmosphere for downward-going 
neutrinos, and altitude of the r decay in the case of Earth-skimming v T - 
spans a sufficiently wide range of numerical values as to guarantee that 
at the edges of the phase space none of the simulated showers fulfills the 
global trigger conditions. This is taken as a clear indication that a complete 
sample of showers has been produced without introducing any bias and 
therefore that the Monte Carlo sample correctly represents the characteristic 
of showers that could trigger the SD of the Pierre Auger Observatory. For the 
Earth-skimming channel, showers were simulated at zenith angles between 
90.1° and 95.9° and at an altitude of the decay point above the Pierre 
Auger Observatory up to 2500 m. In the case of downward-going neutrinos, 
simulations were performed at zenith angles in the range 75°-89°. 

5. Identifying neutrino-induced showers 

As stated above, the selection of potential neutrino-induced showers 
(neutrino candidates) is based on two steps: 

1. Firstly, we select among the data collected at the SD of the Pierre 
Auger Observatory those events that arrive in inclined directions with 
respect to the vertical. 
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Figure 4: Sketch of the different types of showers induced by UHE neutri- 
nos. All the channels depicted contribute to the neutrino event rate due to 
downward-going v induced showers. 

2. Secondly, we select among the inclined events those with FADC traces 
that are spread in time, indicative of the presence of an inclined shower 
in the early stage of development, a clear signature of a deeply inter- 
acting neutrino triggering the SD. 

Although the two steps above are the same for all the neutrino-induced 
showers searched for at the Pierre Auger Observatory, due to the different 
nature of Earth-skimming and downward-going neutrino induced showers, 
the criteria and selection cuts that are applied to data are slightly different. 



5.1. Selection of inclined events 

First of all, events occurring during periods of data acquisition instabil- 
ities [25] are excluded. 

For the remaining events the FADC traces of the triggered stations are 
first "cleaned" to remove accidental signals induced (mainly) by atmospheric 
muons arriving closely before or after the shower front - produced in showers 
different than the triggering one and which are below the energy threshold 
of the Pierre Auger Observatory. The trace cleaning procedure is detailed 
in |27j . After that, the start times of the signals in all stations included 
in the global trigger are requested to be compatible with a plane shower 
front moving at roughly the speed of light. This compatibility is realized 
through upper bounds on both, the largest residual and the mean quadratic 
residual from the planar fit. If the condition is not fulfilled, fits are attempted 
removing one station; for this operation, the stations are sorted by increasing 
quality (based on the integrated amplitude and the duration of the signal) , 
and the procedure is stopped as soon as a satisfactory solution is found. 
If none is found, trials are made removing two stations, and so on. The 
event is accepted if at least three (four) stations in the Earth-skimming 
(downward-going) case belong to the configuration. 

The second step in both channels is the selection of inclined showers. 
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From the pattern (footprint) of stations at ground (see Fig.[5|) we can extract 
a length L along the arrival direction of the event (i.e., the main axis of the 
event) and a width W perpendicular to it characterizing the shape of the 
footprint (see [27J for complete details). The ratio L/W depends on zenith 
angle. Vertical events have L/W ~ 1 and this ratio increases gradually as 
the zenith angle increases. Very inclined events typically have elongated 
patterns on the ground along the direction of arrival, and hence large values 
of L/W. A cut in L/W is therefore a good selector of inclined events. The 
exact value of this cut is different for downward-going and Earth-skimming 
events, and was determined through Monte Carlo simulations of ^-induced 
showers performed at different zenith angles. For downward-going events 
with 8 > 75° the requirement is L/W > 3, while for Earth-skimming it is 
more restrictive L/W > 5 since only quasi-horizontal showers with largely 
elongated footprints can trigger the array u (see Fig. 3 in Ref.|27]). 




Figure 5: Schematic view of the footprint of a shower triggering the surface 
detector array of the Pierre Auger Observatory. The shower triggers the ar- 
ray from the left to the right of the figure, along the "main axis" . The circles 
represent the position of the stations, with their sizes being proportional to 
the collected signal in the PMTs. See text for more details. 

Another indication of inclined events is given by the apparent speed 
V of the trigger from a station j to a station j, averaged over all pairs 
(i,j) of stations in the event. This observable denoted as (V) is obtained 
in a straightforward manner from the distance between the stations after 
projection along the "main axis" of the footprint at ground (dij) as depicted 
in Fig. [5l and from the difference in trigger times of the stations (Atij). 
Vertical showers have apparent average speeds exceeding the speed of light 
since all triggers occur at roughly the same time, while in very inclined 



1 The axis of Earth-skimming showers travelling in the upward direction does not 
intersect ground, contrary to the downward-going showers case. For this reason, we exploit 
the properties of the footprint generated by the shower particles that deviate laterally from 
the shower axis and trigger the water Cherenkov stations. 
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events (V) is concentrated around the speed of light. Moreover its Root- 
Mean-Square (RMS(V)) is small. For downward-going (Earth-skimming) 
events (V) is required to be below 0.313 m ns _1 ((V)G [0.29,0.31] m ns _1 ) 
and RMS(y)/(V}< 0.08 (RMS(y)< 0.08 mns' 1 ). The values of these 
selection requirements are based on comparisons between data and Monte 
Carlo simulations. Also, and only for downward-going events, a further 
quality cut is applied consisting on a simple reconstruction of the zenith 
angle rec and the requirement that 6> rec > 75° (see [28] for full details). 

In the top part of Table Q] the cuts applied to the observables used to 
select inclined events are summarized. 





Earth-skimming 


Downward-going 




Number of Stations > 3 


Number of Stations > 4 


Inclined 
Showers 


L/W >5 

0.29 m ns" 1 <(V)< 0.31 m ns" 1 

RMS(F) < 0.08 m ns" 1 


#rcc > 75° 

L/W >3 
(V) < 0.313 m ns- 1 
RMS(V)/(V)< 0.08 


Young 
Showers 


At least 60% of stations with 
ToT trigger & AoP > 1.4 


Fisher discriminant T based 
on Area-over- Peak (AoP) 



Table 1: Observables and numerical values of cuts applied to select inclined 
and young showers for Earth-skimming and downward-going neutrinos. See 
text for explanation. 



5.2. Selection of young showers 

Once inclined showers are selected the next step is to identify young 
showers among the data collected at the SD of the Pierre Auger Observatory. 

To optimize the numerical values of the cuts and tune the algorithms 
needed to separate neutrino-induced showers from the much larger back- 
ground of hadronic showers, we divided the whole data sample into two 
parts (excluding periods of array instability). A fraction of the data (train- 
ing period) is dedicated to define the selection algorithm. These data are 
assumed to be overwhelmingly constituted of background showers. The 
applied procedure is conservative because the presence of neutrinos would 
result in a more severe definition of the selection criteria. The remaining 
fraction is not used until the selection procedure is established, and then it 
is "unblinded" to search for neutrino candidates. In Table [2] we indicate the 
periods used for training and "blind" search. The blind search period for the 
Earth-skimming (downward-going) analysis corresponds to an equivalent of 
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Earth-skimming 


Downward-going 


Training period 
Blind search period 


1 Nov 04 - 31 Dec 04 
1 Jan 04 - 31 May 10 


1 Jan 04 - 31 Oct 07 

1 Nov 07 - 31 May 10 


Equivalent full Auger 
blind search period 


3.5 yr 


2.0 yr 


v candidates 








Diffuse limit 90% C.L. 
(GeV cm -2 s _1 sr _1 ) 


k < 3.2 x 10~ 8 


k < 1.7 x 10~ 7 


Energy range (EeV) 


~ 0.16-20.0 


~ 0.1 - 100.0 



Table 2: Training and blind search periods for the search for Earth-skimming 
and downward-going neutrino candidates. In the 3 rd row we indicate the 
equivalent period of time of a full surface detector array. In the 4 row we 
give the number of candidates found in the search period after unblindly 
applying the cuts selecting inclined and young showers (see Table [1]) . In the 
5 th row we give the numerical value of the 90% C.L. limit to the normal- 
ization A; of a diffuse flux of UHE neutrinos assumed to behave with energy 
as dN/dE = k E~ 2 . Systematic uncertainties are included in the value of 
the limit (see Sec. 16.31 for details). In the last row we indicate the energy 
range where the limits apply, typically the energy interval where 90% of the 
events are expected. 



~ 3.5 yr (~ 2 yr) of a full surface detector array consisting of 1600 stations 
working continuously without interruptions. 

It is worth remarking that data instead of Monte Carlo simulations of 
hadronic showers are used to optimize the identification cuts. The first 
reason for this is that, the composition of the primary UHECR flux - a 
necessary input in the simulations - is not accurately known. Also, the 
detector simulation may not account for all possible detector defects and/or 
fluctuations that may induce events that constitute a background to UHE 
neutrinos, while they are accounted for in collected data, including those 
which are not well known, or even not yet diagnosed. 

This is the general strategy followed in the search for Earth-skimming 
v T and downward-going z^-induced showers. However the two searches differ 
in several aspects that we detail in the following sections. 

5.2.1. Earth-skimming analysis 

In the Earth-skimming analysis we identify young showers by placing 
a cut on the fraction of stations in the event that fulfill two conditions: 
(1) the station passes the ToT local trigger condition, and (2) the ratio of 
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the integrated signal over the peak height - the so-called Area-over-Peak 
(AoP), a variable that carries information on the time spread of the signal 

- is greater than 1.4. By convention, both the "area" and the "peak" values 
are normalized to 1 in signals induced by isolated muons. 

The aim of both conditions is to identify broad signals in time such 
as those induced by showers developing close to the array. In particular, 
with the second condition we reject background signals induced by inclined 
hadronic showers, in which the muons and their electromagnetic products 
are concentrated within a short time interval, exhibiting AoP values close 
to the one measured in signals of isolated muons. 

In order to reject inclined hadronic events, at least 60% of the triggered 
stations in the event are required to fulfill the two conditions above (Table[T]). 
The selection conditions were optimized using data collected during the 
training period indicated in Table It is important to remark that this is 
the same selection procedure and training period as in previous publications 
[26} 127] . which is applied in this work to a larger data set. The final choice 
of the actual values of the neutrino selection cuts was done by requiring zero 
background events in the training data sample. When the Earth-skimming 
cuts in Table Q] are applied blindly to the data collected during the search 
period, no events survived. 

5.2.2. Downward-going analysis 

In the search for downward-going events, the discrimination power is 
optimized with the aid of a multi-variate technique known as the Fisher 
discriminant method [34J. The method consists on constructing a linear 
combination of observables denoted as J- which optimizes the separation 
between two samples of events, in our case background hadronic inclined 
showers occuring during the downward-going training period (see Table [2]) , 
and Monte Carlo simulated v- induced showers. The method requires as 
input a set of variables which can discriminate between the two samples. 
For that purpose we use variables depending on the Area-over-Peak (AoP) 

- as defined above - of the FADC traces. In the first few stations hit by a 
deep inclined shower, the typical AoP values range between 3 and 5 (left- 
hand panel of Fig. [6]). 

After training the Fisher method, a good discrimination is found when 
the following ten variables are used to construct the linear Fisher discrim- 
inant variable T: the AoP of the four stations that trigger first (early sta- 
tions) in each event, their squares, the product of the four AoPs, and a 
global parameter that measures the asymmetry between the average AoP of 
the early stations and those triggering last (late stations) of the event. 

The product of the AoP of the earliest four stations in the event aims 
at minimizing the relative weight of an accidentally large AoP produced, 
for instance, by a single muon which does not belong to the shower front 
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Figure 6: Distributions of the Area-over- Peak (AoP, see text) of the earliest 
station (left) and of the product of the AoP of the first four stations in the 
event (right). In each panel we show the distribution of the corresponding 
variable in background events (i.e. data events in the training sample as 
indicated in Table [2]) , and in simulated electron neutrino charged-current 
events. These are two of the ten variables depending on the AoP used 
in constructing the multi-variate Fisher discriminant linear polynomial to 
optimize the separation between background and neutrino-induced showers. 
See text for more details on the remaining 8 variables. 

arriving at a station before or after the shower itself. This variable is also a 
very good discriminator as shown in the right-hand panel of Fig. [6l We have 
also checked in Monte Carlo simulations that neutrino-induced events typ- 
ically have an asymmetry parameter larger than proton or nucleus-induced 
showers. 

As the shower front is broader at larger distance from the core for both 
young and old showers, the discrimination is better when splitting the sam- 
ples according to the number of selected stations N. A Fisher discriminant 
polynomial was obtained separately for 4 < N < 6, 7 < N < 11, and 
N > 12. An excellent separation is achieved for events in each of the three 
sub-samples. The individual AoPs of the first four tanks have the largest 
weights in the Fisher polynomials. In Fig. [7] we show as an example the 
distribution of J- in the sub-sample with the smallest number of selected 
stations (the distributions corresponding to the three sub-samples can be 
found in Fig. 7 of [28J). 

Once the Fisher discriminant T is defined, the next step is to define a 
numerical value of T, denoted as T cu t, that separates neutrino candidates 
from regular hadronic showers. One of the advantages of the Fisher discrimi- 
nant method is that it allows us to estimate the expected rate of background 
events, and hence to tune the value of J-" cu t so that the background is kept 
at a very low value. This is important given the fact that the expected 
rate of detected neutrino events will be small. Data in the training period 
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Figure 7: Distribution of the value of the Fisher polynomial (J-, see text 
for details) for events with number of selected stations 4 < N < 6. Data in 
the training period (see Table [2D describe the nucleonic background, while 
Monte Carlo simulated downward-going neutrinos correspond to the signal. 
The vertical lines indicate J- C ut needed to expect 1 event in the labeled 
periods of time (full SD array) . 



indicated in Table [2] was exploited to produce a reasonable prediction of the 
background (see [28J for full details). In practice, we fix J- cn t so that the 
estimated number of background events is 1 in 20 yr of data taking by a full 
Auger SD. With this cut, and for our search sample we have an estimated 
background of 0.1 events for each multiplicity class that add up to a total 
of 0.3 events with a statistical uncertainty of 30%. It is important to re- 
mark that this estimate relies on the a priori hypothesis that the background 
has an exponential distribution in T . Given the fact that we do not have 
a solid estimation of the actual background, a conservative approach was 
taken assuming the background is zero, in other words, the estimated 0.3 
background events were not used to improve our upper limit on the flux [37J 
(see Sec. O below). 

As exemplified in Fig. [7] for the low multiplicity events, the identifica- 
tion cuts reject only ~ 10% of the simulated neutrino events, and those 
are mainly neutrinos interacting far from the ground that, being similar to 
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nucleonic-induced showers, are not expected to be identified. 

Finally, when the downward-going cuts in Table Q] are applied to the 
data collected during the search period, no neutrino candidates appeared 
(see Table E|). 

6. Exposure to UHE neutrinos 

6.1. Neutrino identification efficiencies 

With the criteria to select neutrino-induced showers indicated in Table [H 
we obtain a relatively large identification efficiency both for Earth-skimming 
v T and downward-going i/-induced showers. The efficiency has been com- 
puted with Monte Carlo simulations as the fraction of simulated events 
identified as neutrinos. 

In the case of Earth-skimming v T induced showers, and a full Auger SD 
working without interruption, the efficiencies depend only on the energy 
of the emerging r leptons (E T ) and on the altitude of the "center of the 
shower" (h c ) above ground (averaged over the decay channels). This is 
conveniently defined as the altitude of the shower axis at a distance of 10 
km away from the r decay point along the shower axis. Showers induced 
by r leptons with the same energy but with different zenith angles - the 
range in 9 being very narrow - have approximately the same efficiency as 
long as the corresponding altitudes of their shower maxima h c are the same. 
The maximum efficiency that can be reached is 82.6%, the 17.4% remaining 
corresponds to the channel in which the r decays into a /i which is unlikely to 
produce a detectable shower close to ground. In Fig. [8] we show the trigger 
and identification efficiencies as a function of h c for different r energies. 
As expected, the efficiency increases with E T and drops as the r decays at 
increasing altitude from ground. 

In the case of downward-going neutrinos the identification efficiency de- 
pends on neutrino flavour, type of interaction (CC or NC), neutrino energy 
(E u ), zenith angle (9), and distance (D) measured from ground along the 
shower axis at which the neutrino is forced to interact in the simulations. An 
example of the efficiency that can be achieved in a full SD array is shown in 
Fig. El The efficiency is different from zero between a minimal depth close 
to ground (a minimal amount of matter needed for the ^-induced shower 
to reach a sufficient lateral expansion), and a maximal one (such that the 
electromagnetic component is almost extinguished at ground level and hence 
the neutrino cannot be identified). The efficiency as well as the slice of at- 
mosphere where it is different from zero, typically increase with neutrino 
energy, and depend on the neutrino flavour and interaction. As an extreme 
example, high energy v T interacting in the atmosphere through the CC chan- 
nel can be identified regardless of the interaction depth in the atmosphere, 
as long as the energetic r produced in the interaction decays and produces 
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a shower close to ground. 

6.2. Exposure 

Ideally, for the calculation of the exposure of the SD of the Auger Obser- 
vatory to ultrahigh energy neutrinos, the simulated neutrino showers should 
be randomly distributed over the actual configurations of the array, applying 
to the shower at ground the trigger and neutrino identification conditions 
to obtain the active (effective) area of the array at every second, and as a 
function of the parameters of the neutrino-induced showers (neutrino energy, 
zenith angle, h c ,...). A sum over time and integration in solid angle would 
then yield the exposure {£) to UHE neutrinos in both the Earth-skimming 
and downward-going neutrino analyses. During the search periods consid- 
ered for both Earth-skimming and downward- going neutrino searches, the 
surface detector array of the Pierre Auger Observatory was growing contin- 
uously. Since the number of working stations and their status are monitored 
every second, we know with very good accuracy the SD configuration at any 
instant as well as its evolution with time. 

In practice, to avoid having to cope with an unaffordable number of 
configurations, different strategies were devised to calculate in an accurate 
and less time-consuming manner the effective area of the SD array to Earth- 
skimming and downward-going z^-induced showers. 

For downward-going neutrinos, the calculation of the exposure involves 
folding the SD array aperture with the v interaction probability and the 
identification efficiency, and integrating in time. Changes in the configura- 
tion of the array introduce a dependence of the efficiency e on the position 
of the core of the shower r = (x, y) in the surface S covered by the array, 
and on time t. 

Assuming a 1:1:1 flavour ratio (as expected due to the effects of neutrino 
oscillations during propagation from the sources), the total exposure can be 
written as [28] : 

S DG (E V ) = —Y a\E v ) [ dtdddD 

i L J 

smOcosO Aig(9,D,E v ,t)] (1) 

where the sum runs over the three neutrino flavours and the CC and NC 
interactions, with a % the corresponding ^-nucleon interaction cross-section 
[35j and m the nucleon mass. The integral is performed over the zenith 
angle 6, the interaction depth D of the neutrino (in units of g cm" 2 ), and 
the blind search period. A l eS is the effective area of the SD array given by: 

Ai s (E u , 6, D,t) = J e\r, 6, D, E v , t) dA (2) 
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where the integral is performed over the core positions r of the showers. 

For the Earth-skimming neutrinos the calculation of the exposure is de- 
scribed in [27] . 

The exposures obtained for the search periods indicated in Table [2] are 
plotted in Fig. [TUJ where for the downward-going neutrino induced showers 
we also plot the contribution of the different channels (Fig. H]) to the total 
exposure. Among them we have included the possibility that downward- 
going v T interact with the mountains surrounding the Observatory which 
provide a dense target for neutrino interactions. 

The exposure to Earth-skimming neutrinos is higher than that to down- 
ward-going neutrinos by a factor between ~ 2 and ~ 7 depending on the 
neutrino energy, partially due to the longer search period in the Earth- 
skimming analysis ~ 3.5 yr of full Auger, compared to ~ 2.0 yr in the case 
of the downward-going analysis. When normalized to the same search time, 
the Earth-skimming channel is still a factor ~ 2.5 — 3 more sensitive when in- 
tegrated over the whole energy range, mainly due to the larger density of the 
Earth's crust where v T interactions can occur, compared to the atmosphere. 
The larger number of neutrino flavours and interaction channels that can 
be identified in the downward-going analysis, as well as the broader angular 
range (75° < 9 < 90° compared to 90° < 8 < 95°), partly compensate the 
difference. 

6.3. Systematic Uncertainties 

Several sources of systematic uncertainty have been carefully considered. 
Some of them are directly related to the Monte Carlo simulation of the 
showers, i.e., generator of the neutrino interaction either in the Earth or in 
the atmosphere, parton distribution function, air shower development, and 
hadronic model. Others have to do with the limitations on the theoretical 
models estimating, for instance, the interaction cross-section or the r energy 
loss at high energies. Some of these sources play a dominant role on the 
Earth-skimming analysis, while others do on the downward-going neutrino 
one. 

In both analyses the procedure to incorporate the systematic uncertain- 
ties is the same. Different combinations of the various sources of systematic 
uncertainty render different values of the exposure, and the final uncertainty 
is incorporated in the value of the limit itself through a semi-Bayesian ex- 
tension [36] of the Feldman-Cousins approach [37]. In Table [3] we summa- 
rize the dominant sources of systematic uncertainty and their impact on 
the exposure. In the Earth-skimming analysis the model of energy loss for 
the r is the dominant source of uncertainty, since it determines the energy 
of the emerging rs after propagation in the Earth; the impact of this on 
the downward-going analysis is much smaller since r energy losses are only 
relevant for u T interacting in the mountains, a channel that is estimated to 
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Source of uncertainty 


Earth-skimming 


Downward-going 


Monte Carlo simulation of shower 


+20%, -5% 


+9%, -33% 


z^-nucleon cross-section 


+5%, -9% 


+7%, -7% 


r energy losses 


+25%, -10% 


+6%, -6% 


Topography 


+18%, 0% 


- 



Table 3: Main sources of systematic uncertainty and their impact on the 
Earth-skimming [27] and downward-going [28] exposures. 

contribute only ~ 15% to the total exposure. The uncertainty on the shower 
simulation, that stems mainly from the different shower propagation codes 
and hadronic interaction models that can be used to model the high energy 
collisions in the shower, contributes significantly in both cases. The presence 
of mountains around the Observatory - which would increase the target for 
neutrino interactions in both cases - is explicitly simulated and accounted for 
when obtaining the exposure of the SD to downward-going neutrino-induced 
showers, and as a consequence does not contribute directly to the systematic 
uncertainties. However, it is not accounted for in the Earth-skimming limit 
shown in Table [21 Instead, we take the topography around the Observatory 
as a source of systematic uncertainty and we estimated that accounting for 
it would have increased the event rate by ~ 18% (Table [3]). 

7. Results 

We have searched for neutrino candidates over the search data periods 
and no events fulfilling either the Earth-skimming or the downward-going 
selection cuts were found. This allows us to put limits to the UHE diffuse 
neutrino flux. 



7.1. Limits to the diffuse flux of UHE neutrinos 

Under the assumption that the UHE neutrino flux <&(E) behaves with 
neutrino energy as: 



.(*)-§ -nr 



[GeV^ 1 cm" 2 s" 1 sr" 1 ] 
the integrated limit on the value of k is: 



N, 



up 



(3) 



(4) 



/£™ x E~ 2 £{E) dE 

where £(E) is the exposure. The actual value of the upper limit on the signal 
events (-/V up ) depends on the number of observed and expected background 
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events. We recall here that, according to [37], N up = 2.44 at 90% C.L. 
for zero candidates and no expected background events. When systematic 
uncertainties are included (Sec. 16. 3p the value of N np changes. 

The final limits are reported in Table [2] where we give the normalization 
k obtained in the search periods (indicated in the same table) for the Earth- 
skimming and downward-going searches. 

In Fig. [TT]we show the Earth-skimming and downward-going integrated 
neutrino flux which indicate the level of a diffuse neutrino flux assumed to 
behave with energy as E~ 2 , needed to detect N up events with a Poisson 
probability of ~ 90% given the exposure accumulated during the 3.5 years 
for Earth-skimming neutrinos (2.0 years for downward-going) of equivalent 
time of a full SD. 

Another way of presenting the results is to display the upper limit in 
differential form. In this procedure we assume that the diffuse neutrino flux 
behaves as E~ 2 within energy bins of 0.5 width on a decimal logarithmic 
scale, and is given by 2.44/(0.5 log(10) • E ■ 6(E)), assuming again no back- 
ground. The differential limit obtained in this way is shown in Fig. [11] for 
the Earth-skimming and downward-going cases. We achieve most (~ 90%) 
of the sensitivity in the energy range ~ 0.16 — 20 EeV (~ 0.1 — 100 EeV) 
for Earth-skimming (downward-going) neutrinos. In Fig. [TTJ we also show 
several predictions of different theoretical models of cosmogenic neutrino 
production [30j[TO]. Predictions for cosmogenic neutrino fluxes depend on 
several unknown parameters including the evolution with redshift of the 
sources and the injected UHECR composition. Given the uncertainties in 
these parameters, and in particular the possible presence of heavy primaries 
in the UHECR spectrum [51], we have plotted a range of models to illustrate 
the wide range of predictions available [TO] . 



7.2. Event rate predictions 

In Table H] we give the expected number of events from a diffuse flux 
of cosmogenic neutrinos (produced in the interaction of cosmic ray protons 
with background radiation fields) |40| . from a model of neutrino production 
through the bottom-up mechanism in Active Galactic Nuclei (AGN) [42] . 
and from a theoretical model |43] in which neutrinos are the product of the 
decay of super-heavy relic particles of the early stages of the Universe. Opti- 
mistic theoretical flux predictions for cosmogenic neutrinos are within reach 
of our present sensitivity and some models of neutrinos produced in accel- 
erating sources are already being constrained. Exotic models are severely 
disfavored. Note that all such top down models are also tightly constrained 
by the limits of the Pierre Auger Observatory on the photon fraction in 
UHECR 
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Model & reference 


Earth-skimming 


Downward-going 


Cosmogenic (Fermi) [40J 

AGN 02] 

Exotic (SH relics) 03] 


- 0.6 

~ 5.1 
~ 3.0 


~ 0.1 
~ 0.8 
~ 1.0 



Table 4: Number of expected events for several theoretical models of UHE 
neutrino production, given the exposure of the surface detector array of the 
Pierre Auger Observatory to Earth-skimming and downward-going neutri- 
nos (Table [2D- 

8. Summary and Prospects 

In this paper we have reviewed the searches for astrophysical sources of 
ultrahigh energy neutrinos at the Pierre Auger Observatory [261 12"?] |2"B] . 

The neutrino detection technique is based on the observation of exten- 
sive air showers induced by downward-going neutrinos of all flavours as they 
interact with the atmosphere, and by upward-going z/ T 's through the Earth- 
skimming mechanism. These ^-induced showers display characteristic fea- 
tures that allow us their identification in the overwhelming background of 
regular UHE hadronic showers. At ground level, high zenith angle neu- 
trino events would have a significant electromagnetic component leading to 
a broad time structure of detected signals in the surface detector array, in 
contrast to nucleonic-induced showers. 

We have shown that, using Monte Carlo simulations and training data 
samples, identification criteria for UHE neutrinos can be defined and used 
to perform a blind search on the remaining data sample. The analysis of the 
collected data at the Pierre Auger Observatory until 31 May 2010 reveals no 
candidate events for either downward- going or Earth-skimming neutrinos. 
Based on this negative result, stringent limits have been placed on the diffuse 
flux of UHE neutrinos. Even though the Auger Observatory was designed 
to measure properties of UHECRs, the limits reported in Table [2] provide at 
present one of the most sensitive bounds on neutrinos at EeV energies, which 
is the most relevant energy to explore the predicted fluxes of cosmogenic 
neutrinos. 

There are several lines of work in progress inside the Auger Collaboration 
related to the neutrino search which will be the subject of future reports. 
Some of the efforts concentrate on the combination of the downward-going 
and Earth-skimming channels into a single analysis. This will simplify the 
search procedure and will obviously translate into an improvement of the 
diffuse neutrino limit. The extension of the downward-going neutrino search 
to lower zenith angles {6 < 75°) is also very promising. Exploring the sky 
down to 6 ~ 60° implies a sizeable increase on the exposure and hence on 
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the limit in case no candidates are found. The main drawback of decreasing 
9 is that the atmosphere slant depth reduces and nucleonic-induced showers 
look "younger" when arriving at ground, making their separation from u- 
induced showers more challenging. On the other hand, the sensitivity to 
neutrino detection could also be extended to lower energies by reducing the 
separation between SD stations. Monte Carlo studies indicate that using a 
configuration of stations similar to the currently existing "infill" array (~ 60 
stations spaced by 750 m) would lead to a significant increase of the neutrino 
detection probability at lower energies (below 0.3 EeV) with respect to the 
standard SD array. Nevertheless, due to the small size of the current infill 
array, the exposure does not appear to be competitive. 

Finally, it is worth mentioning that the sensitivity of the Pierre Auger 
Observatory to the detection of UHEj/s from potential astrophysical point- 
like sources is being evaluated. The absence of candidates in the searches 
for diffuse neutrino fluxes described in this report allows us to place limits 
on the neutrino fluxes coming from sources in the field of view of the SD of 
the Auger Observatory. Preliminary results indicate that with the SD we 
are sensitive to a large fraction of the sky spanning ~ 100° in declination 
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Figure 8: T3 trigger (open dots) and identification (closed dots, cuts as 
in Table [1]) efficiency in the Earth-skimming analysis, as a function of the 
height above ground of the shower at 10 km from the r decay point h c . 
The efficiency is shown for Monte Carlo showers induced by rs with energy 
(clockwise from the top left panel) 0.1, 1, 10 and 100 EeV. The efficiencies 
are calculated in a full SD array (see text for details). 
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Trigger and identification efficiency for v e CC channel: 85° - 1 EeV 
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Figure 9: Fraction of electron neutrinos of energy 1 EeV and 6 = 85° trig- 
gering the array (solid grey line) and passing the downward-going analysis 
cuts in the second column of Table Q] (solid black line) as a function of the 
slant depth of the interaction above the ground. The dashed line represents 
the fraction of events passing all cuts except for the cut on the Fisher dis- 
criminant J- (see Sec. 15. 2|) . The efficiencies are calculated in a full SD array 
(see text for details). 
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Figure 10: Exposure of the surface detector array of the Pierre Auger 
Observatory on the data search periods to Earth-skimming z^-induced show- 
ers (equivalent to 3.5 yr of full Auger), and to downward-going z/-induced 
showers (equivalent to 2 yr of full Auger) . 
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Figure 11: Thick lines: Differential and integrated upper limits (at yuyo 
C.L.) to the diffuse flux of UHE neutrinos (single flavour assuming equipar- 
tition) from the Pierre Auger Observatory for downward-going v (equivalent 
search period = 2 yr of full Auger) and Earth-skimming v T (equivalent search 
period = 3.5 yr of full Auger). Limits from other experiments are also plot- 
ted |38l 139] . All limits have been scaled to single flavour. The IceCube 
differential limit is scaled by a factor 1/2 due to the different binning in 
energy with respect to the Auger differential limits. 

Thin lines: Expected fluxes for three theoretical models of cosmogenic neu- 
trinos (scaled to single flavour when necessary), "p, Fermi-LAT" [10] corre- 
sponds to the best fit to UHECR spectrum incorporating Fermi-LAT bound 
assuming that the transition from Galactic to extragalactic CRs takes place 
at 10 19 eV. "p, evol-FRII" [10] assumes the FRII strong source evolution 
with a pure proton composition, dip transition model and maximum energy 
of UHECRs at the sources E Ptmax = 10 2L5 eV. "Fe, uniform" [TO] represents 



an extreme model assuming an iron rich composition, low E pmax , 
evolution of the UHECR sources case. 



uniform 
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